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Abstract 

A  ceiling-mounted  semi-cylindrical  model  containing  a  single  wrap-around  fin 
(WAF)  was  tested  in  the  AFIT  Mach  2.9  test  facility.  Flow  visualization  using  oil- 
fiow  streaklines,  schlieren  images  and  shadowgraph  photography  revealed  a  A-shock 
at  the  fin-body  juncture  and  the  development  of  an  asymmetric  bow-shock  about  the 
fin.  Quantitative  measurements  were  taken  with  a  10°  cone-static  pressure  probe,  a 
Pitot  pressure  probe  and  two  cross-wire  hot-film  probes  (u-v  and  u-w  components, 
respectively).  Measurements  were  made  at  cutting-planes  from  the  inlet  of  the  test 
section  to  aft  of  the  model,  with  emphasis  placed  in  the  vicinity  of  the  WAF.  Results 
include  cutting-plane  profiles  and  contours  of  mean  and  turbulent  fluctuations  of 
the  primitive  and  conserved  flow  variables.  It  was  found  that  the  incompressible 
turbulent  fluctuating  quantities  are  equally  as  descriptive  of  the  flow  structure  in 
the  fin’s  vicinity  as  the  compressible  turbulence  fluctuations.  The  asymmetric  bow- 
shock  was  found  to  be  an  inviscid  phenomenon  which  was  stronger  on  the  concave 
side  than  the  convex  side  and  of  diminishing  strength  at  the  tip  with  no  bleeding 
effects  over  the  tip. 
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MACH  2.9  INVESTIGATION  INTO  THE  FLOW 
STRUCTURE  IN  THE  VICINITY  OF 
A  WRAP-AROUND  FIN 


I.  Introduction 

1.1  Problem  Statement 

Current  military  design  emphasis  has  been  placed  on  low  observable  properties 
for  advanced  fighter  aircraft.  To  reduce  the  intensity  of  radar  return  signatures 
of  airframes,  researchers  are  investigating  the  use  of  internal  weapon  stowage  and 
decreasing  the  use  of  sharp  corners.  A  design  solution  to  these  problems  is  to  utilize 
wrap-around  fins  (WAFs)  on  the  missiles  carried  by  these  aircraft.  A  WAF  is  a  lift 
or  control  surface  which  has  the  same  curvature  of  the  cylinder  to  which  it  is  mated. 
The  key  feature  of  a  WAF  is  that  when  stowed  it  fits  flush  with  the  cylinder  surface. 

Wrap-around  finned  missiles  (with  stowed  fins)  can  be  stored  in  a  compartment 
of  less  volume  than  a  missile  body  of  the  same  size  which  utilizes  straight,  non¬ 
folding  fins.  With  WAF  missiles  employed,  an  airframe  can  carry  more  missiles  in 
the  volume  previously  used  by  the  straight  fin  missiles.  Alternatively,  the  airframe 
can  use  the  internal  volume  savings  for  fuel  or  avionics.  Wrap-around  fin  missiles 
simplify  the  design  of  the  airframe  housing  and  launching  equipment  on  airframes 
which  integrate  the  weapon  in  partial  submersion.  Partial  submersion  WAF  designs 
avoid  difficulties  associated  with  fin-fuselage  contact.  An  additional  feature  of  the 
WAF  missile  is  that  the  fin  stowed  configuration  has  a  lower  radar  return  due  to 
the  loss  of  the  fin-body  juncture.  Therefore,  a  low  observables  benefit  is  achieved 
by  using  WAF  missiles  on  external  pylons  for  current  operational  aircraft.  While 
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missile  and  airframe  designers  are  interested  in  WAFs  for  the  previously  mentioned 
benefits,  there  are  interesting  features  of  WAF  missile  designs. 

There  are  three  key  differences  between  WAF  missiles  and  straight-finned  mu¬ 
nitions:  WAFs  experience  a  rolling  moment  reversal  near  Mach  1,  a  coupling  between 
pitch  and  yaw,  and  a  possible  dynamic  instability  associated  with  undesirable  side 
moments.  These  WAF  performance  characteristics  are  coupled  to  the  flow  structure 
about  the  fin.  At  the  time  of  this  publication,  a  detailed  characterization  of  the 
flowfield  in  the  vicinity  of  a  WAF  has  not  been  made.  This  study  represents  an 
attempt  to  experimentally  address  this  void  in  the  literature. 

A  natural  first  step  in  looking  at  the  WAF  flowfield  is  to  eliminate  the  influence 
of  other  fins.  The  present  effort  considers  a  single  WAF  attached  to  a  partial  cylinder, 
based  on  the  dimensions  of  the  USAF  basic  research  missile  (BRM)  model.  The  fin 
is  placed  at  0°  angle-of  attack  to  create  a  baseline  case.  However,  the  “real  world” 
missile  effects  are  not  entirely  removed.  The  WAF  is  mounted  to  half  of  a  semi¬ 
cylinder,  to  place  the  fin  out  of  the  natural  boundary  layer  of  the  tunnel  while 
providing  for  the  investigation  of  the  curvature  effects  associated  with  missile-fin 
junctures. 

An  additional  goal  of  this  research  is  to  provide  information  of  value  to  compu¬ 
tational  fluid  dynamics  (CFD)  research.  To  that  end,  conventional  hot-wire  anemom- 
etry  measurements  provide  the  necessary  mean  flow  and  turbulence  information  re¬ 
quired  to  validate  numerical  efforts.  Great  care  is  taken  in  detailing  the  upstream 
conditions  of  the  flowfield,  to  enable  numerical  duplication  and  expansion  of  this 
work.  Additionally,  the  influence  of  turbulence  and  compressibility  issues  on  numer¬ 
ical  simulation  is  examined.  It  is  the  intent  of  the  author  that  this  information  will 
help  guide  further  experimental  and  numerical  research  in  this  area. 
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Fin  Detail 

Figure  1.1  USAF  Basic  WAF  Research  Model  (WL-TR-94-7015) 

1.2  Background 

WAFs  find  uses  in  many  military  applications  -  from  tube  launched  munitions 
for  the  Army  and  Navy,  to  internally  stowed  weapons  in  the  Air  Force.  These 
multiple  applications  produce  an  array  of  WAF  designs  to  study.  The  WAF  chosen 
for  this  design  is  one  patterned  after  the  BRM  (7).  Figure  1.1  shows  the  pattern 
of  the  BRM  with  a  wrap-around  fin.  The  BRM  differs  from  the  WAF  model  by  a 
straight  fin  with  the  same  planform  area. 

Numerous  experimental  tests  are  accomplished  in  instrumented  ballistic  ranges 
(75,  41,  10,  12,  71,  11,  68).  The  majority  of  these  tests  cover  the  subsonic  and 
transonic  flight  regime  (75,  41,  10,  12,  11).  These  tests  show  a  roll  moment  reversal 
near  Mach  1.0,  a  coupling  between  yaw  and  pitch  and  a  potential  for  side  moments 
leading  to  undesirable  dynamic  instabilities.  The  supersonic  experiments  tend  to 
show  a  loss  of  static  stability  near  Mach  4.0  (71).  Although  these  range  experiments 
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provide  valuable  design  information,  the  facilities  are  not  designed  for  detailing  the 
flow  structure  of  the  WAF. 

Some  wind  tunnel  tests  on  WAFs  have  been  conducted  (8,  47).  In  the  paper  by 
Kretzschmar,  et  al.  (47),  the  focus  was  on  designing  lift  and  control  surfaces  and  not 
of  flowfield  description.  Abate  and  Hathaway  (8)  provided  very  useful  information 
on  the  surface  structure  of  the  missile  with  multiple  WAFs.  Lacking  from  the  sum  of 
these  previous  experiments  was  information  pertaining  to  the  characteristics  of  the 
shock  structure  about  the  WAF  and  turbulence  information  near  the  WAF. 

Numerical  simulations  were  performed  on  missile  configurations  with  a  view 
toward  characterizing  the  shock  structure  (9,  25,  54,  70).  The  study  of  the  Hyper¬ 
sonic  Applied  Research  Technology  (HART)  missile  showed  a  pressure  bleeding  effect 
over  its  straight  fins  at  angle-of- attack  (54).  The  remaining  studies  were  performed 
on  wrap-around  fins  and  did  not  model  turbulence. 

Settles  and  Dodson  (61)  have  defined  criteria  necessary  in  experimental  data 
gathering  for  CFD  research  in  supersonic  turbulent  shock-boundary  layer  interaction. 
A  characteristic  of  this  flowfield  associated  with  the  WAF  configuration  tested  is  the 
shock/boundary-layer  interaction  near  the  juncture  of  two  curved  surfaces.  There¬ 
fore,  these  criteria  provide  appropriate  guidance  for  the  current  research  effort.  Using 
their  criteria.  Settles  and  Dodson  have  identified  a  number  of  high-quality  experi¬ 
mental  datasets  suitable  for  CFD  research  into  turbulence  model  validation.  Each 
of  these  experiments  are  supersonic  “building-block”  experiments  containing  turbu¬ 
lence  information.  Some  studies  in  the  database  have  considered  the  interaction 
between  straight  fins  and  flat  plates  (43,  37,  49).  The  current  study  represents  an 
incremental  advance  over  typical  shock  generation  schemes  by  introducing  surface 
curvature  effects.  The  primary  emphasis  of  creating  these  criteria  is  the  need  to  nu¬ 
merically  simulate  the  experimental  environment.  To  that  end,  this  effort  contains 
a  rigorous  treatment  of  the  conditions  at  the  inlet  to  the  test  section.  Additionally, 
the  uncertainties  associated  with  the  current  study  are  defined  and  quantified. 
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1.3  Objectives 

The  goals  of  this  experimental  study  are  two-fold.  The  primary  goal  is  to  ex¬ 
perimentally  characterize  the  mean  and  turbulent  flowfield  in  the  vicinity  of  a  wrap¬ 
around  fin.  The  shock  structure  is  determined,  using  a  myriad  of  flow  visualization 
techniques;  shadowgraph  imagery,  schlieren  photography  and  oil-flow  surface  streak¬ 
lines.  Additionally,  mean  flow  quantities  and  turbulence  information  are  measured 
by  Pitot,  cone  static  pressure  and  cross- wire  hot-film  anemometers  in  the  flowfield 
about  the  WAF.  An  alternate  goal  is  to  provide  this  information  in  a  useful  form  for 
validation  of  numerical  turbulence  modeling. 

1.4  Outline 

This  section  provides  a  road  map  of  the  chapters  which  follow.  In  Chapter 
II,  the  facility  for  this  experiment  was  detailed  along  with  the  methodology  of  hot¬ 
wire  measurements,  shadowgraph  and  schlieren  photography,  and  oil-flow  streakline 
visualization.  The  process  used  in  the  reduction  of  measured  voltages  to  valuable 
information  about  the  primitive  and  conserved  flow  variables  is  detailed  Chapter  III. 
The  results  of  the  flowfield  survey  are  presented  in  Chapter  IV.  Finally,  Chapter  V 
provides  a  series  of  the  conclusions  and  suggests  the  direction  for  future  efforts. 

Additional  information  is  included  in  the  appendices.  Appendix  A  details 
the  uncertainty  in  measurements  and  the  effect  of  these  uncertainties  on  the  re¬ 
sults.  Appendix  B  provides  the  coordinates  of  the  AFIT  Mach  2.9  Test  Facility’s 
converging- diverging  nozzle. 
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II.  Experimental  Apparatus  and  Methodology 

This  chapter  details  the  equipment  and  techniques  used  to  perform  the  exper¬ 
imentation  of  this  thesis. 

2.1  Mach  2.9  Facility 

The  AFIT  Mach  2.9  test  facility  is  an  intermittent  high  pressure  blow-down 
tunnel  with  downstream  evacuation.  The  tunnel  is  located  in  Building  640,  Room 
148  on  Wright-Patterson  AFB,  Area  B.  This  facility  features  a  run  time  of  30  seconds 
with  an  evacuation  time  of  6  to  10  minutes.  Each  test  section  measures  6.35  x  6.35  x 
33.02cm.  Figure  2.1  is  a  photo  of  the  nozzle  and  the  two  removable  test  sections. 


Figure  2.1  Nozzle  and  Test  Sections  with  WAF  Model 


The  coordinate  system  applied  throughout  this  thesis  was  the  hot-film  coordi¬ 
nate  system.  The  origin  for  this  system  was  defined  to  be  at  the  nozzle  throat  on  the 
floor’s  port  side  (looking  downstream).  The  x-axis  runs  longitudinally  downstream 
and  the  y-axis  is  positive  up.  The  z-axis  runs  horizontally  towards  the  starboard  wall 
(facing  downstream)  -  thus  completing  the  right-handed  orthogonal  triad.  Figure  2.2 
locates  the  coordinate  axes  at  the  origin. 

The  vacuum  tanks  (with  a  volume  of  ~16m^)  are  evacuated  by  one  Stokes 
Microvac  Pump  -  model  212-11,  and  two  IPP  Evacuation  Vacuum  Pumps  -  model 
FW-EV-PI  each  powered  by  a  Reliance  XE  Dutymaster  A/C  motor  requiring  5.6kW. 
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The  minimum  vacuum  achievable  is  ~2mmHg,  though  a  typical  run  required  only 
lOmmHg.  Figure  2.3  shows  the  vacuum  pump  system. 

The  high  pressure  air  is  delivered  at  6.8atm  from  the  house  air  supply,  which  is 
powered  by  two  Atlas  Copco  GAU  807  compressors  delivering  0.45^.  The  moisture 
is  removed  from  the  air  by  two  Pioneer  R500A  refrigeration-type  air  dryers.  The  air 
flows  to  the  tunnel  from  the  compressors  and  dryers  by  ~150m  of  7.62cm  diameter 
cast  iron  pipe.  Due  to  a  pressure  drop  during  the  course  of  a  tunnel  run,  the  pressure 
is  regulated  by  a  Leslie  Regulator  Valve  class  GPK-1,  size  2  -  which  is  controlled 
by  the  Fairchild  Pressure  Regulator  System  model  HP  190A.  The  air  passes  through 
the  pressure  regulator  into  the  plenum  chamber  and  flow  straightener,  which  are 
shown  in  Figures  2.4  and  2.5  -  respectively.  The  plenum  and  flow  straightener  are 
connected  to  the  nozzle  in  Figure  2.6. 

The  air  flow  is  accelerated  through  a  two-dimensional  converging- diverging 
nozzle  as  described  in  Figure  2.7.  The  coordinates  originate  from  the  throat  on  the 
tunnel  floor.  Nozzle  coordinates  are  tabulated  in  Appendix  B. 

At  this  point,  the  flow  enters  sequentially  through  two  removable  test  sections 
and  into  the  vacuum  tanks  by  way  of  an  adjustable  diffuser.  Figure  2.8  is  a  schematic 
of  the  facility. 
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Figure  2.3  Vacuum  Pump  System 


Figure  2.4  Plenum  Chamber 
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Figure  2.6  Tunnel  Layout  -  Upstream  Section 
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Figure  2.7  Mach  2.9  Nozzle 
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Figure  2.8  Schematic  of  AFIT  Mach  2.9  Test  Facility  with  WAF  Model 


At  the  entrance  to  the  test  section,  the  freestream  Mach  number  is  2.85  with  a 
standard  deviation  of  1.27%  and  a  measurement  uncertainty  of  1.8%.  Figure  2.9  is  a 
contour  plot  of  the  Mach  number  at  the  inlet  to  the  test  section.  Note  the  Cartesian 
coordinates  are  non-dimensionalized  by  the  fin  radius. 


test  section  ceiling 


Qi  ■  '  '  '  I  I  ■  ■  ■  i  .  ■  .  .  I  I  ■  ■  ■  I 

^012,34 

z/r 


Figure  2.9  Mach  Contours  at  the  Test  Section  Inlet 

The  boundary  layer  thickness  is  5.3mm  (based  on  95%  of  Moo)  at  a  distance 
31.1cm  from  the  throat  -  just  downstream  of  the  test  section  entrance  -  with  an 
uncertainty  of  4.1%  in  measured  position.  The  settling  chamber  operates  nomi¬ 
nally  with  a  total  pressure  and  total  temperature  of  2.07atm  and  297K,  yielding  a 
freestream  Reynolds  number  of  1.70  x  W /m.  The  freestream  turbulence  (defined 
as  the  square  root  of  two-thirds  the  non-dimensional  compressible  turbulent  kinetic 
energy  (TKE)  (74))  is  0.8%  with  a  standard  deviation  of  0.2%.  The  compressible 
turbulent  kinetic  energy  at  the  test  section  entrance  is  shown  in  Figure  2.10. 

The  uncertainties  related  to  the  documentation  of  measured  freestream  quan¬ 
tities  are  explained  in  Appendix  A. 
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Figure  2.10  Compressible  TKE  Contours  at  the  Test  Section  Inlet 


2.2  Wrap-Around  Fin  Model 

The  WAF  model  is  a  cylinder  of  the  fin  radius  blended  to  the  removable  test 
section  ceiling  at  a  peak  height  of  0.79cm.  The  cylinder  is  8.14cm  in  length  with 
the  single  fin  placed  at  the  downstream  base.  Upstream,  the  cylinder  is  blended 
to  the  tunnel  floor  with  a  polynomial  spline  which  ensures  second  order  continu¬ 
ity.  The  blending  function  (Figure  2.11)  is  translated  into  physical  coordinates  in 
Equation  2.1. 


C  —  677®  —  15?/^  +  IO77® 

7]  =:  0.2f  -  4.296,  21.48  <  f  <  26.48 
^1.0 -(^-2.0)2 -0.5,  1.134  <  ^  <  2.866 

K 


^  =  4.0 


(2.1) 


The  blending  region  is  7.94cm  long  and  starts  34.10cm  from  the  nozzle  throat.  The 
fin  has  a  thickness  of  0.32cm,  a  span  of  2.25cm  and  a  chord  length  of  2.03cm.  The 
leading  edge  and  tip  of  the  fin  are  beveled  at  45°.  Figure  2.12  is  a  perspective  view 
of  the  WAF  model  in  the  test  section  ceiling.  Figure  2.13  is  an  orthographic  drawing 
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with  dimensions  -  a  detail  of  the  WAF  model.  The  model  is  then  inserted  into  a  test 
section  ceiling  as  shown  in  Figure  2.14. 


Figure  2.12  WAF  Model 

The  floor  of  the  test  section  was  designed  with  a  series  of  interchangeable  plugs. 
The  plugs  were  5.08cm  wide  and  0.64cm  thick  with  the  ends  rounded  by  a  0.32cm 
radius.  These  plugs  allowed  the  0.64cm  diameter  probe  sleeve  to  move  across  the 
floor  in  the  z-axis.  A  large  plug  measuring  5.08  x  3.81cm  contained  a  slot  2.54cm 
long  and  0.64cm  wide  -  with  all  corners  rounded  by  a  0.32cm  radius.  The  probe 


Figure  2.14  WAF  Model  in  Test  Section  Ceiling 
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sleeve  can  slide  across  the  y-axis  near  the  fin  (during  operation  the  sleeve  is  secured 
by  a  pressure  plate).  Figure  2.15  is  a  photo  of  the  test  section  floor  and  removable 
plugs. 


Figure  2.15  Test  Section  Floor 


2.3  Flow  Visualization 

Shadowgraphs  and  schlieren  images  were  taken  using  Xenon  Corporation  model 
437B  nanopulser,  with  a  spark  duration  of  10  nanoseconds,  as  a  light  source  (6).  The 
images  were  recorded  on  Polaroid  type  57  high  speed  instant  black  and  white  film. 
For  shadowgraphs,  the  light  source  was  aimed  away  from  the  test  section  to  be 
reflected  back  by  a  1.524m  focal  length  mirror  -  placed  1.524m  from  the  center  of 
the  test  section.  It  can  be  shown  that  the  intensity  of  a  shadowgraph  correlates 
to  the  second  order  gradient  of  index  of  refraction  in  the  plane  normal  to  the  light 
beam  (32).  The  index  of  refraction  is  related  linearly  to  the  density  of  the  medium  the 
light  passes  through.  Figure  2.16  is  a  schematic  of  the  shadowgraph  setup.  Photos 
were  taken  through  both  the  Plexiglas  and  the  optical  grade  glass.  The  optical  grade 
shadowgraphs  were  of  much  higher  quality  than  the  Plexiglas  photos.  Figure  2.17 
shows  the  optical  grade  glass  walls  installed  on  the  first  test  section. 
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Camera  Body 


Figure  2.16  Shadowgraph  Schematic 
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The  schlieren  setup  introduces  another  concave  mirror,  with  1.534m  focal 
length  -  as  well  as  a  knife  edge.  The  light  passing  through  the  test  section  from 
the  shadowgraph  is  now  focused  and  half  of  the  image  is  blocked  by  the  knife  edge 
at  the  focal  point.  The  image  is  then  reflected  by  a  flat  mirror  and  focused  through  a 
lens  -  to  cope  with  space  limitations.  The  intensity  of  the  recorded  schlieren  image  is 
equivalent  to  the  gradient  of  density  normal  to  the  knife  edge  (32).  Therefore,  a  knife 
edge  positioned  horizontally  will  display  the  vertical  density  gradient.  Rotating  the 
knife  edge  180°  will  produce  the  negative  image.  Figure  2.18  is  the  schlieren  setup. 


Camera  Body 

Figure  2.18  Schlieren  Schematic 

Both  photographic  systems  were  aligned  using  an  air-bubble  torpedo  level. 
First  the  test  section  was  leveled.  Then  the  mirrors  were  bore-sighted  using  an  Oriel 
He-Ne  gas  laser  model  79255  producing  a  0.48mm  diameter  linearly-polarized  beam 
with  a  wavelength  of  632. 8nm  and  a  power  of  0.5mW  (3).  No  focusing  was  needed 
for  the  shadowgraph.  However,  the  schlieren  was  focused  by  using  the  light  source  to 
read  an  overhead  transparency  on  a  ground-glass  plate  inserted  in  the  camera  body. 
The  transparency  was  suspended  in  the  center  of  the  tunnel  test  section  through  one 
of  the  removable  plugs. 
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Surface  streakline  patterns  were  visualized  using  an  oil  mixture  of  stamp  ink 
and  Dow  Corning  200  Fluid,  1000  cSt.  General  purpose  cleaning  fluid  was  used  to 
mix  the  ink  and  oil.  For  these  runs,  the  WAF  model  was  placed  in  the  test  section 
floor  -  to  keep  the  oil  from  dripping  off  before  the  test.  The  oil  mixture  was  applied  to 
the  cylinder,  floor  and  fin.  Then  the  tunnel  was  run  for  15  to  20  seconds.  The  results 
were  photographed  after  the  tunnel  had  stopped  -  allowing  for  the  removal  of  the 
walls.  The  tunnel  unstart  caused  negligible  changes  to  the  observed  flow  patterns. 

8.4  Pressure  and  Temperature  Measurements 

All  pressure  measurements  were  made  by  Endevco  pressure  transducers.  These 
transducers  measure  gauge  pressure.  The  ambient  pressure  was  measured  at  the  be¬ 
ginning  and  end  of  a  series  of  runs  -  about  every  two  hours  -  on  a  Wallace  &  Tiernan 
model  233111  absolute  pressure  gauge.  The  Pitot  pressure  in  the  test  section  was 
measured  with  a  Pitot  probe.  Figure  2.19.  The  cone-static  pressure  was  measured 
with  a  cone-static  probe  which  had  a  10°  ±  0.03°  semi-vertex  angle  as  seen  in  Fig¬ 
ure  2.20.  The  probe  pressure  was  measured  by  an  Endevco  model  8510C-15.  The 
signal  was  processed  through  an  Endevco  model  4423  signal  conditioner  and  pow¬ 
ered  by  an  Endevco  model  4225  power  supply.  The  voltage  output  was  displayed 
on  an  HP  3446A  digital  multimeter.  Just  downstream  of  the  flow  straightener,  the 
plenum  pressure  is  measured  through  a  Pitot  tube  with  an  Endevco  pressure  trans¬ 
ducer,  model  8510C-100  and  is  reported  by  an  Endevco  pressure  display  unit,  model 
4428A.  The  plenum  temperature  was  measured  by  an  Omega  K  type  thermocouple 
and  displayed  on  an  Omega  model  15B  display  unit  (5). 

2.5  Hot-Film  Anemometry 

Constant  temperature  cross- wire  hot-film  anemometry  measurements  were  used 
in  the  majority  of  the  flowfleld  analysis.  The  employed  system  consisted  of  the  TSI 
IFAIOO  Intelligent  Flow  Analyzer  which  contained  three  model  150  constant  tem- 
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Figure  2.20  Cone-Static  Pressure  Probe 
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perature  anemometers  with  a  series  resistance  (Rg)  of  50f2.  Two  anemometers  were 
used,  each  attached  to  one  channel  of  the  cross-wire  probe.  The  two  cross-wire  probes 
used  were  the  model  1243-20  (u-v  plane)  and  1243AN-20  (u-w  plane)  hot-film  probes. 
Each  probe  has  a  similar  configuration,  differing  only  by  the  orientation  of  the  wires. 
The  wires  were  actually  made  of  a  thin  film  of  platinum  51/um  in  diameter  and  1mm 
long.  The  measured  resistances  of  the  ‘wire’  and  attached  cable  are  summarized  in 
Table  2.1 


Table  2.1  Measured  Wire  Resistances, 


Probe  Serial  Number 

Channel  1  (0) 

Channel  2  (0) 

934011 

5.898 

5.822 

934012 

6.088 

6.184 

944011 

6.880 

6.887 

934012 

6.010 

6.898 

Overheat  ratios  of  ~2.0  were  set  by  an  external  breadboard  and  high  quality 
resistors,  in  combinations  of  If)  and  lOf).  Great  care  was  taken  in  placement  of 
the  resistors  on  the  breadboard  in  an  attempt  to  reduce  inductance  and  resulting 
oscillatory  behavior.  The  cross-wire  probes  were  seated  in  a  TSI  two-sensor  probe 
support  for  placement  into  the  tunnel.  Figures  2.21  and  2.22  are  schematics  of  the 
u-v  and  u-w  hot-film  probes,  respectively. 


Figure  2.21  U-V  Cross-Wire  Hot-Film  Probe  Schematic  (TSI,  1987) 

To  properly  determine  the  flowfield  properties,  the  hot-film  response  must  be 
calibrated.  The  calibration  was  accomplished  by  placing  the  probe  in  the  test  section 
at  a  point  of  known  Mach  number  and  varying  the  Reynolds  number.  The  calibration 
point  for  this  study  was  at  the  center  of  the  test  section  and  31.12cm  downstream 
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Figure  2.22  U-W  Cross-Wire  Hot-Film  Probe  Schematic  (TSI,  1987) 

of  the  nozzle  throat.  The  Mach  number  used  in  calibration  was  2.9.  (The  measured 
Mach  number  at  the  calibration  point  was  2.83,  a  2.4%  difference.)  The  Reynolds 
number  was  varied  by  adjustments  to  the  plenum  pressure.  A  curve  fit  was  obtained 
by  least  squares  regression  relating  the  wire  voltage  from  the  anemometer  to  the 
local  Reynolds  number.  A  more  detailed  mathematical  explanation  of  the  calibration 
technique  is  found  in  Chapter  III. 

2.6  Data  Acquisition 

Quantitative  measurements  were  made  at  165  individual  stations  throughout 
the  test  section.  These  stations  formed  y-z  cutting  planes.  Each  station  was  referred 
to  by  two  reference  positions.  The  first  reference  was  the  non-dimensional  position 
The  second,  and  more  interesting,  reference  was  to  the  non-dimensional  position 
The  latter  reference  position  was  the  distance  from  the  nozzle  throat  to  the  plane 
in  which  the  unflexed  hot-film  probe  is  located.  Note  that  the  Pitot  and  cone-static 
probes  were  actually  longer  than  the  hot-film  probe,  by  6.35mm.  The  difference  in 
lengths  resulted  in  a  determination  of  Mach  number  which  was  located  upstream  of 
the  hot-film  probe  measurements  made  from  the  same  station.  This  difference  in 
position  was  accounted  for  in  the  region  near  the  fin,  where  the  pressure  probes  were 
aligned  to  the  hot-film  probe.  However,  the  measurements  made  upstream  of  the 
fin  implicitly  assumed  a  small  change  in  pressure  over  the  distance  between  pressure 
probe  and  hot-film  probe  tips.  Figure  2.23  shows  the  y-z  cutting  planes  formed  by 
the  hot-film  tip,  which  include  flex  effects.  The  ^  reference  position  is  found  by 
extrapolating  the  position  in  Figure  2.23  down  to  the  floor.  The  reference  positions 
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(^)  associated  with  Pitot,  cone-static  and  hot-film  probe  tips  corresponding  to  each 
cutting  plane  in  Figure  2.23  are  detailed  in  Table  2.2. 


Figure  2.23  Hot-Film,  Y-Z  Cutting  Plane  Locations 


Table  2.2  Reference  Positions 


Reference 

Letter 

Unflexed  Cross- Wire 
Probe  Position 
(Reference) 

Unflexed  Pressure 
Probe  Position 

X 

X 

r 

r 

A 

19.6 

19.2 

B 

26.8 

26.4 

C 

29.8 

29.4 

D 

30.0 

29.6 

E 

30.2 

29.8 

F 

30.4 

30.0 

G 

30.6 

30.2 

H 

30.8 

30.4 

I 

31.0 

30.6 

J 

31.2 

30.8 

K 

31.4 

31.0 

L 

31.6 

31.2 

The  flowfield  was  mapped  by  attaching  a  particular  probe  to  a  custom  made 
traverse  stand  pictured  in  Figure  2.25.  The  stand  has  a  Unislide  Assemblies  series 
A2500  traverse  slide  with  a  thread  pitch  of  40  per  2.54cm  (Figure  2.24).  To  digitally 
record  the  traverse  position,  a  Linear  Displacement  Voltage  Transducer  (LDVT) 
was  attached  between  the  moving  slide  and  the  traverse  stand.  The  LDVT  was  pow- 
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ered  from  a  Hewlett  Packard  6205C  Dual  DC  Power  Supply.  The  output  from  the 
LDVT  was  run  through  the  LDVT  Signal  Conditioner  model  1000-0012  to  produce 
the  recorded  voltage.  The  traverse  slide  assembly  can  be  adjusted  manually,  or  by 
computer  with  the  help  of  a  stepper  motor  and  Arrick  Robotics  model  MD-2  Dual 
Stepper  Motor  Driver.  The  driver  was  controlled  by  Microsoft’s  Quick  BASIC  soft¬ 
ware  and  subroutines  provided  from  Arrick  Robotics  on  a  Compaq  Deskpro  80386DX 
running  at  25MHz. 

SERIES  A2500 

\ 


1  ^  dia 


Figure  2.24  Traverse  Slide  (Unislide  Inc.,  Units  in  Inches) 


Before  starting  the  tunnel,  the  traverse  was  positioned  to  put  the  probe  0.254cm 
away  from  the  ceiling  of  the  tunnel.  This  small  gap  provided  a  measure  of  safety  to 
keep  the  probes  from  brushing  the  ceiling  during  tunnel  operation,  as  the  probe  flexes 
back  as  the  flow  rushes  past.  Probe  position  was  determined  by  ‘sighting’  through 
two  matching  grid  patterns  -  one  affixed  to  each  Plexiglas  tunnel  wall.  As  the  tunnel 
was  run,  the  flexing  of  the  probe  was  very  noticeable.  Measurements  were  made  to 
determine  the  flexing  of  the  probe,  the  results  from  which  are  presented  in  Table  2.3. 
Also,  since  the  cross-wire  probes  determine  velocity  information,  the  flex  angle  of 
the  u-v  was  required  (the  flex  of  the  u-w  probe  had  an  insignificant  effect  on  the 
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angle  each  wire  made  to  the  flow).  The  flex  angle  of  the  u-v  probe  was  determined 
to  be  2.86°. 


Table  2.3  Measured  Probe  Flexing 


Probe 

Ai/ 

y 

Ax 

V 

Pitot 

1.00% 

7.25% 

Cone- Static 

1.00% 

5.00% 

Cross- Wire 

2.25% 

6.25% 

The  general  strategy  applied  to  mapping  the  flowfleld  was  to  insert  a  probe 
into  one  of  the  slots  in  the  tunnel  floor  and  take  measurements  with  the  traverse  in 
motion.  The  traverse  speed  was  0.19^.  The  freestream  flow  speed  at  the  entrance  to 
the  test  section  was  608.1  j.  Therefore,  the  traverse  speed  introduced  only  negligible 
errors  to  the  data  acquisition.  Taking  data  while  the  traverse  was  in  motion  allowed 
almost  the  entire  test  section,  4.5cm,  to  be  covered  in  one  tunnel  run  -  greatly 
speeding  the  acquisition  process.  At  the  end  of  one  tunnel  run,  the  probe  would 
be  positioned  0.32cm  away  from  the  previous  position  within  the  slot  on  the  tunnel 
floor.  This  process  was  repeated  until  the  entire  plane  was  covered.  The  mapping  of 
every  plane  was  accomplished  with  one  probe,  before  moving  onto  the  next  probe. 
For  example,  the  Pitot  probe  was  used  to  map  the  tunnel  before  proceeding  onto 
the  cone-static  probe  and  cross-wire  probes.  This  strategy  was  altered  slightly  for 
hot-fllm  mapping.  The  calibration  was  done  at  the  beginning  or  end  of  a  mapping 
plane.  Since  the  probe  had  to  be  removed  to  exchange  slots  in  the  tunnel  floor,  the 
probe  was  placed  at  the  calibration  point  for  one  calibration  run  and  then  placed 
in  the  next  slot  in  the  sequence.  The  calibration  procedure  was  determined  to  be 
acceptable  given  the  variation  in  calibration  coefficients  between  runs  was  negligible 
as  seen  in  Table  2.4. 

The  plenum  pressure,  probe  pressure,  or  hot-film  voltage  and  position  were 
recorded  on  a  Nicolet  MultiPro  digital  acquisition  (D/A)  unit  containing  four  model 
120  digital  acquisition  boards  -  each  capable  of  sampling  four  channels  at  IMHz. 
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i  Calibration  Constants 


mnel  1 

Channel  2 

a 

b 

Overheat 

Ratio 

a 

b  ' 

0.178 

mtumm 

2.01 

0.172 

-0.360 

0.173 

2.01 

0.166 

-0.303 

liWWl 

-0.397 

|I2|£| 

-0.316 

0.181 

-0.432 

0.176 

-0.270 

0.181 

-0.379 

0.183 

-0.290 

1.98 

0.182 

-0.406 

-0.330 

1.98 

0.183 

-0.431 

0.183 

-0.323 

1.98 

-0.409 

2.00 

0.149 

2.00 

0.164 

-0.378 

0.151 

-0.362 

2.00 

0.158 

0.156 

-0.367 

2.00 

0.160 

-0.389 

2.00 

0.162 

IilEEl 

-0.394 

2.00 

0.166 

-0.406 

torn 

-0.428 

2.00 

lilWl 

-0.409 

bm 

-0.413 

2.00 

liMlifel 

-0.377 
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The  D/A  unit  was  controlled  by  Microsoft  Windows  software  running  on  a  Zenith 
80486DX  33MHz  personal  computer.  The  plenum  pressure,  Pitot  and  cone-static 
pressure  and  position  from  the  LDVT  were  sampled  at  200Hz,  while  the  anemometer 
voltages  were  sampled  at  lOkHz.  A  faster  sampling  rate  was  desired  for  the  hot-film 
probes;  however,  the  D/A  unit  would  only  store  256Kbytes  of  data.  At  lOkHz, 
26.2  seconds  of  data  could  be  stored.  The  traverse  required  23.4  seconds  to  span 
the  tunnel.  Therefore,  a  faster  sampling  rate  would  require  multiple  runs  to  span 
one  line  of  data.  159  individual  lines  of  data  were  mapped  in  this  experiment, 
requiring  almost  1300  tunnel  runs  and  250  hours  of  tunnel  time  -  which  constrained 
the  sampling  rate  to  lOkHz.  Since  the  Nyquist  frequency  of  this  sampling  rate  is 
5kHz,  no  spectral  decomposition  was  performed  and  the  data  was  not  background 
filtered  at  60Hz. 
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III.  Data  Reduction 


Though  the  Navier-Stokes  equations  govern  the  motion  of  Newtonian  fluids, 
modifications  have  been  made  to  the  Navier-Stokes  equations  to  account  for  the 
wide  range  of  length  and  time  scales  associated  with  turbulent  flowfields.  In  such, 
approximate  forms  of  the  primitive  variables,  such  as  velocity,  density,  pressure  and 
temperature,  can  be  split  into  two  terms 


Q  —  Q  -\-  Q' 


Typically,  the  first  term  is  a  time  averaged  quantity  and  the  second  term  is  the  fluc¬ 
tuation  about  the  mean  value.  From  this  variable  separation,  the  Reynolds- averaged 
Navier-Stokes  (RANS)  equations  are  formed.  However,  the  mean  and  fluctuating 
values  obtained  from  hot-wire  anemometry  do  not  directly  correlate  to  the  variables 
used  in  RANS  (with  the  exception  of  Reynolds  turbulent  shear  stress).  Significant 
manipulations  (the  scope  of  this  chapter)  are  required  to  produce  the  normalized 
fluctuations  common  to  RANS  analysis. 

CFD  applications  typically  involve  the  use  of  conservative  variables,  such  as 
mass-flux  and  energy.  CFD  tasks  use  mass-weighted  time  averaging,  which  is  related 
to  RANS  variables  by 

Q  =  Q  +  Q" 

Q  =  f_ 

qTi  ^ 

These  relationships  are  the  basis  of  the  Favre-averaged  Navier-Stokes  (FANS)  equa¬ 
tions.  It  is  important  to  note  that  hot-wire  results  must  also  be  manipulated  to 
provide  data  of  value  for  CFD.  The  separation  of  variables  into  RANS  and  FANS 
form  requires  the  use  of  controversial  assumptions  and  the  omission  of  higher  order 
terms.  These  simplifying  assumptions  degrade  the  accuracy  of  the  results.  Nonethe- 
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less,  these  results  are  intended  to  give  CFD  analysts  a  comparison  for  many  of  the 
variables  which  they  generate. 

This  chapter  examines  the  method  used  to  reduce  data  obtained  via  the  Pitot 
pressure  probe,  cone-static  pressure  probe,  u-v  and  u-w  hot-film  cross-wire  probes. 
Mean  flow  information  such  as  density,  pressure,  and  velocity  magnitude  are  derived 
from  the  pressure  probes.  The  hot-film  cross-wire  probes  determine  velocity  vectors 
and  turbulence  values  in  a  form  which  are  manipulated  to  be  useful  for  comparative 
analysis  to  CFD. 

3.1  Data  Sample  Reduction 

Up  to  25  Megabytes  of  data  was  collected  per  run  of  the  wind  tunnel.  Run¬ 
ning  the  wind  tunnel  1300  times  produced  a  prohibitively  large  amount  of  data  for 
analysis.  To  avoid  wasting  scarce  computer  resources,  the  relevant  information  from 
these  data  files  was  collected  and  stored.  The  technique  used  to  extract  the  relevant 
information  follows.  Mean  values  of  large  samples  were  used  to  reduce  both  pressure 
and  hot-film  records.  The  mean  is  determined  by 

1  ^ 

^  ~  lu  (^-1) 

where  N  is  the  sample  size  and  x  is  the  recorded  information.  Additionally,  the 
fluctuation  information  from  hot-film  samples  was  found  with  the  following  relations 

(3.2) 

_  1  ^ 

^'y'  ^  jfZTi  -  y)  (3.3) 

where  the  former  equation  is  the  variance  and  the  latter  equation  is  the  covariance 
between  two  variables,  x  and  y.  These  sample  reduction  techniques  reduced  a  total 
900  Megabytes  of  collected  data  to  11  Megabytes  of  stored  information. 
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3.2  Pressure  Probes 


The  measurement  of  total  pressure  behind  a  normal  shock,  combined  with 
the  measurement  of  static  pressure  behind  a  shock  created  from  a  10°  cone,  allows 
determination  of  the  local  Mach  number  in  a  given  flow.  The  total  pressure  behind 
a  normal  shock  {Pt2)  is  measured  directly,  using  a  Pitot  probe.  The  static  pressure 
behind  a  shock  from  a  10°  cone  (Pcs)  is  measured  using  a  10° ±0.03°  cone-static  probe. 
Through  the  manipulation  of  results  from  (44,  14)  and  normal  shock  relations,  the 
following  curve  fit  was  generated  by  Bowersox  (18) 

^  =  -0.052976  ±  4.6840^  -  18.6786^^  ±  50.7006^^  -  54.1577^^  (3.4) 


where 

>  _  -^ti  ^cs 

^  ~  Pt2  Pn 

Note  that  the  plenum  pressure  (Pn)  is  introduced  to  minimize  the  differences  between 
the  two  separate  tunnel  runs  required  to  get  the  Pitot  and  cone-static  probe  data. 
It  has  been  found  for  flow  angles  less  than  6.0°  that  the  errors  in  Mach  number  are 
less  than  0.3.  Also  note  that  Equation  3.4  is  valid  for  Mach  numbers  in  the  range 
from  1.5  to  4.4  and  has  a  standard  deviation  of  0.06%. 

3.3  Hot-Film  Probes 

This  section  examines  multiple  overheat  (MOH)  hot-wire  data  reduction  of 
Bowersox  (17,  18)  as  an  intermediate  step  to  obtaining  the  single  overheat  (SOH) 
reduction  equations. 

The  constant  temperature  hot-wire  anemometer  records  the  voltage  required  to 
maintain  the  wire  at  a  constant  known  temperature.  The  power  required  to  maintain 
this  temperature  is  equivalent  to  the  heat  transfer,  between  the  hot-wire  and  the 
surrounding  flow.  The  Nusselt  number,  the  non-dimensional  heat  transfer  coefficient. 
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can  be  related  to  the  heat  transfer  from  the  wire  by 


irkL{Tw  -  Te) 

where  k  is  the  coefficient  of  thermal  conductivity  for  the  fluid,  L  is  the  wire  length, 
Tw  is  the  wire  temperature  and  Te  is  the  equilibrium  temperature  or  the  temperature 
the  unheated  wire  would  approach  under  these  specific  flow  conditions.  The  wire 
temperature  is  calculated  from 

Tw  =  Tref  H - ^  {Rw  —  Rref)  (3-5) 

Ter  Ttrej 

where  Tref  is  the  reference  temperature,  is  the  temperature  coefficient  of  resistance 
(0.0024//t')  and  Rref  is  the  reference  wire  resistance  at  Tref  -  The  wire  heat  transfer 
rate  is  equal  to  the  power  supplied  to  the  hot-film,  or 


ijR 


w 


where  is  the  wire  voltage,  i-ei  is  the  current  supplied  to  the  hot-wire  and  R^  is 
the  wire  resistance.  The  wire  current  can  be  related  to  the  wire  voltage  by 

.  gw 

Rw  T  Rs  +  Rl 

where  Rg  and  Rl  are  the  series  and  lead  resistances,  respectively.  The  coefficient  of 
thermal  conductivity  is  determined  from  the  following  curve  fit  (in 

k  =  -5.30377  •  10“^  +  1.02108  •  lO’^T*  -  4.70285  •  10"®r/  -f-  1.30279  ■  lO-^^T^^ 


The  final  problem  with  determining  the  Nusselt  number  is  the  difficulty  of  mea¬ 
suring  the  equilibrium  temperature.  However,  where  the  effective  Reynolds  number 
is  greater  than  20,  the  difference  between  Te  and  T  is  less  than  3.0%.  This  error 


3-4 


can  be  eliminated  if  total  temperature  is  consistently  substituted  for  equilibrium 
temperature  in  both  calibration  and  data  reduction.  With  the  above  manipulations 
in  place,  the  Nusselt  number  can  be  found  by  the  following  equation 


Nu  = 


^  2  p 

Oyj  -itxiJ 

irkL^Tui  —  Rg  Rl)^ 


(3.6) 


3.3.1  Calibration.  It  has  been  determined  experimentally  that  Kings  Law, 
the  functional  relationship  between  Nusselt  number  and  effective  Reynolds  number 
for  incompressible  flow,  is  also  an  acceptable  relationship  for  compressible  flow  (45). 
Kings  Law  takes  the  following  functional  form 


Nu  =  a 


rRCg  b 


(3.7) 


The  effective  Reynolds  number  is  defined  in  terms  of  the  cosine  law  by  Spangen- 
burg  (66)  as 

R  A, 

Ree  = - cos<p 

1^ 

where  p  is  the  freestream  density,  U  is  the  magnitude  of  the  velocity,  d  is  the  hot- 
wire’s  diameter,  and  ^  is  the  sweep  angle  of  the  hot-wire  with  respect  to  the  flow. 
To  calibrate  the  hot-film  probe,  the  probe  is  placed  under  known  conditions  and 
the  mass  flow  is  changed  by  varying  total  pressure.  Mass  flow  is  related  to  known 
conditions  by  a  manipulation  of  isentropic  flow  and  normal  shock  relations  to  yield 

pU  =  Pt2M 

where  Pt2  is  the  total  pressure  behind  a  normal  shock  (Pitot  pressure),  M  is  the 
Mach  number,  R  is  the  gas  constant  for  air  and  7  is  the  ratio  of  specific  heats  for 


^  (1  -t-  - — -M^ 

RTi  V  2 


(7  +  1)M2 


7  +  1 

2'yM'^  —  7 
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air.  Viscosity  is  determined  from  Sutherland’s  Law 


1.458  •  IQ-^T/ 
Tt  +  110.40 


(3.8) 


The  remaining  unknowns  in  Equation  3.6  are  a  and  6,  which  are  determined  by  linear 
least  squares  regression. 


3.3.2  General  Theory.  Combining  the  relationship  for  Nusselt  number 
(Equation  3.6),  Kings  Law  (Equation  3.7)  and  substituting  for  viscosity  and  ther¬ 
mal  conductivity  with  their  turbulent  counterparts  yields  the  thermal  anemometry 
response  equation 


(3.9) 


where  the  following  relationships  are  used 


kt=  ko{^y'  ,nfc  =  0.89  ,r,  =  600K  ,  A;,  =  45.6  •  10-3^ 
^«  =  A<o  (^)  ^  ,n^  =  0.77  , //o  =  30.3  •  10"6Pas 


^  _  {Rw  Rs  +  RlY  _jj 

Oq  —  'KIjKq 

Xt-i/j 


Re, 


Eo 


pUd 


COS(f> 


Replace  5^,,  Reg^,  and  7)  by  their  mean  and  fluctuating  components,  using  the 
Binomial  Theorem,  and  retain  only  the  first  order  terms,  to  obtain 


«\/Ree<,  +  b  (T^-Tt) 


(3.10) 
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Solving  for  ^  the  thermal  anemometry  fluctuation  equation  is  given  by 

C-ty 


=  =  /(^) 
\  R-^Co  / 


(3.11) 


where  the  hot-wire  sensitivities  for  constant  temperature  anemometry  are 


/=Mi+-4 

^  \  ayRe 


and  g 


^  (3.12) 


2(T„  -  T,)  2 


Squaring  Equation  3.11  and  averaging  yields 


fi 


2 


+  ‘^fi9i 


Ree/T/X 
RCeo  7f/ 


(3.13) 


Three  overheat  ratios  are  required  to  solve  for  the  three  turbulence  terms.  The 
covariance  equation  is  obtained  by  multiplying  Equation  3.11  by  itself  with  each 
equation  representing  a  separate  wire  response.  Averaging  the  result  provides 


(3.14) 


3.3.3  Coordinate  Transformation.  The  response  of  the  cross- wires  is  inde¬ 
pendent  of  Mach  number  for  M  sin  <f>  >  1.  The  effective  Reynolds  number  must  be 
transformed  to  tunnel  coordinates  through  oblique  shock  relations 


cos  (j> 

s\n<j) 

< 

Re^: 

Ret 

—  sin  (f> 

cos  <f) 

Rcy  ^ 

(3.15) 
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The  effective  Reynolds  number  becomes 


Rse^  —  Re„^  +  kJ^Yiet^ 

—  ^iRe^;^  “b  2j42Re2rRej^  T  ylsRe^^ 


(3.16) 


where  Ai  are  given  by 

Ai  =  cos^  (f>  +  kc  sin^  ^ 

A2  =  {1  —  kc^)  cos  (j)  sin  (j)  (3-17) 

A3  =  k^  cos^  (j)  +  sin^  (f) 


and  kc  is 


kc  =  =fc  (3.18) 

P2 

where  k  =  fcn{Lld).  For  this  experiment,  Ljd  =  19.61,  which  relates  to  k  ^ 
0.50.  However,  since  the  flow  angles  variance  is  small  between  calibration  and  data 
recording  (±10.0°),  k  was  set  to  zero  in  both  calibration  and  reduction.  This  was 
done  consistently  to  calibration  and  data  reduction,  thus  minimizing  the  associated 
error. 


Replacing  Reg^,  Re^,^  and  Rej,^  by  their  mean  and  fluctuating  components, 
applying  the  binomial  theorem,  and  using 


B2  =  t 

B3  =  Ai  ±  2A2R0 


(3.19) 


the  following  can  be  shown 


and 


Ree/A 

RSeo  / 


Bij 


'ReJ 


Re^. 


(3.20) 

j  +  B2j  1 

fRe.;A 

1  RSito  J 

(3.21) 
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where  the  index  j  sums  over  the  two  wires  on  the  cross-wire  probe.  Solving  this  set 
of  equations  and  transforming  into  x  and  y  components 


Jj— 2^5*;^ 


bi  -  Ite.,  " 


„  ■ _ ?eo2  M22 

■4ii/A21  —  A12/A22 


1  R-^eoi  /^ll~Reeo2  /^12 


Re  =  _ 

2Rea;<,  A^ijAw  —  ^22/^12 


Squaring  Equation  3.21  and  averaging  results  in  Equations  3.22  and  3.23.  Note  that 
Di  =  {821! B\\  —  B22IB12)  and  D2  =  {BiifB2i  —  B\2lB22)- 


'Re.;' 


1 


Re 


•Xo 


'Re, 


/\  2 


Vo 


Re 


'Xq 


D2^ 


1  fReJ' 


B2^ 


Re, 


0/1 


1  fReeS 


B21B22  \  Rg, 

2 


Ree;\  /Ree;\  ^  1  (ReJ^ 


eo  /  1  V  I^eeo  /  2 


22 


Re 


eo  /  2} 

(3.22) 


Re.,'\  /Re.;' 


Bii'  I  Re 


■eo  /  1 


B\\B\2  V  Re, 


■eo  /  \  V  Reeo  )  2 


+ 


1  (Re, 


/\2 


eo 


B 


12 


Re 


eo  /  2\ 

(3.23) 


Multiplying  two  forms  of  Equation  3.21  together  -  one  for  each  hot-wire  response 
-  then  averaging  the  resulting  equation,  results  in  the  covariance  response,  Equa¬ 
tion  3.24. 


'R^Re^\  _ 


1 


Rea;„  Rea;„ }  2B11B21 


(3.24) 


To  analyze  the  results  from  the  u-w  hot-film  probe,  replace  v  by  u;,  replace  y 
with  2;  and  change  the  sign  of  the  sine  terms  in  the  transformation  matrix. 


3.3.4  Single  Overheat  Analysis.  Sensitivity  to  total  temperature  fluctua¬ 
tion  is  generally  negligible  for  hot-wires  operating  at  high  overheat  ratios  > 

^  ^ref 

2.0).  The  AFIT  Mach  2.9  test  facility  has  been  found  to  maintain  total  tempera¬ 
ture  fluctuation  below  about  2.0%  (52,  24,  51).  The  following  section  describes  the 
analysis  method  used  for  SOH  measurements. 
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The  Strong  Reynolds  Analogy  is  used  to  correct  for  the  total  temperature  sensi¬ 
tivity  of  the  probe.  When  total-temperature  fluctuations  are  neglected,  Morkovin  (57) 
states  that  the  temperature  fluctuation  can  be  represented  by 


(3.25) 


Bowersox  (17)  postulates  that  when  Pr  ^  I  and  T/  ^  0  then  the  temperature 
fluctuations  can  be  represented  by 


T'  (u'  v'^ 


(3.26) 


where  k  ranges  from  0.3  to  0.5.  Defining  0  =  — k(7  —  1)M^,  where  Rl6  <C  1  and 
■<  1  —  then 

u'  1  \(pu)'  „  fpuVl  ^^27) 


u  1  —  6 


pu 


pu 


^  _  Rq  (pu)'  (pv)' 
u  1  —  9  pu  'pu 

The  total  temperature  fluctuation  has  become  a  function  of  the  mass-flux 


(3.28) 


n  _^pae 
T,  i-e 


pu 


pu 


(3.29) 


Setting  kc  =  D  results  in  an  alteration  of  the  transformation  coefficients  A; 


Ax  =  cos^  (f> 

A2  =  cos  <j)  sin  (j) 
A3  =  sin^  <!> 


(3.30) 


Inserting  the  assumption  that  ^ 


0  alters  Equations  3.13  and  3.14  as  follows 


RejV  1 


(3.31) 
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(3.32) 


V  J  1  V  Ree„  J  2  /1/2  [Su,Ji[^wJ2 

Bowersox,  Miller  and  McCann  (17,  52,  24,  51)  have  determined  experimentally 
that  the  data  reduction  is  unaffected  by  neglecting  the  Strong  Reynolds  Analogy. 
Therefore,  k  is  set  to  1.0  which  implies  P  =  —ad.  The  above  assumptions  are  applied 
and  reduce  the  transformation  to 


2f^Bi2B22  fiBl 

/1/2.R11-B12  /:/2  (-B11-B22  +  B12B21)  f\hB2\B2<2. 


Rej:„'  Re„„' 

Rcio  Reio 


> 


(3.33) 


=  < 


> 


Or  substituting  the  results  from  Equations  3.31  and  3.32  into  Equations  3.22, 
3.23  and  3.24  yields 


(3.34) 
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where  [S'ij]  ^  is 


D'iBiJt 

1  (  Bl  _  Bh 

2B11B21/1"  [bJ^ 

D'tiuft 


Cii/I 


)  1 

( 

^11  ) 

/  2B11B21/2 

^1^12/2 


1 

B2B21B22/1/2 

_ ^ _ r_gll .. , _ ^ 

B11B21/1/2  VBiiBi2C^  B21B22DI 
1 

D'iBiiB22hh 


3.3.5  Separation  of  Turbulence  Variables.  One  of  the  key  features  of  cross¬ 
wire  anemometry  is  that  the  Reynolds  turbulent  shear  stress  can  be  directly  mea¬ 
sured  as  the  negative  of  the  mass-flux  correlation  term  combined  with  the  density 
fluctuation  as  shown  below  (20) 


{puyipv)'  p  (p'Y 

{pu)  \P  J 


(3.35) 


The  density  fluctuation  cannot  be  directly  measured  from  the  cross-wire  probe,  and 
must  be  reduced  by  neglecting  pressure  and  temperature  fluctuations.  Since  this 
term  is  squared  and  multiplied  by  Ro,  it  is  negligible.  Therefore  an  error  in  density 
fluctuation  caused  by  the  assumptions  used  in  data  reduction  will  only  produce  a 
very  small  error  in  Reynolds  turbulent  shear  stress. 


Additional  relationships  for  separation  of  conservative  hot-wire  variables  into 
primitive  turbulence  variables  are 


P 


a  +  l3 


pu 


pu 


Tt\  P' 
Tt  P 


(3.36) 


where 


^  ^  {p^y  _  py_ 

u  p 

^  _  (^  _  o  ^ 

u  pu  p 


«  =  [1  +  ^ 


and  ^  =  (y  —  l)a!Af^ 


(3.37) 

(3.38) 
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Kistler  (42)  suggests  p'  is  proportional  to  u'^.  Therefore,  p'  is  neglected  as  it  is  like 
a  second  order  term.  Note  that  the  pressure  fluctuation  is  multiplied  by  a,  which 
is  near  |  in  the  Mach  regime  of  this  study.  This  results  in  a  small  effect  on  the 
reduction  of  the  density  fluctuation.  Bowersox  and  Schetz  (20)  have  experimentally 
verified  the  validity  of  neglecting  the  effect  of  p'  in  hot-wire  analysis  for  a  Mach  4.0 
free  mixing  layer.  After  neglecting  T/  and  p',  Equation  3.37  reduces  to 

-  =  ^ 

p  a  +  /3  \  “pTiy 

Placing  the  previous  three  equations  in  terms  of  the  already  determined  cross-wire 
conservative  fluctuating  terms,  and  noting  that  1  (neglecting  third  order 

terms)  the  incompressible  fluctuation  equations  become 


fMV 

a  +  ^\\  pu  ) 


(3.40) 
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After  neglecting  total  temperature  and  pressure  fluctuations,  the  cross-correlation 
terms  are  determined  as  follows 
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The  Reynolds  turbulent  shear  stress  is  then  calculated  by 


_  U'v'  p'v'  ^  p'u' 
pv?  'pu  ° 

This  method  of  calculation  is  mathematically  equivalent  to  Equation  3.35.  However, 
it  is  provided  as  a  determination  of  the  errors  generated  in  the  process  of  separating 
the  fluctuating  variables. 


The  variables  conserved  in  the  Favre- averaged  Navier-Stokes  equations  are 
available  for  comparison  as  well.  The  Favre  variables  are  available  by  way  of  the 
following  transformations 


u  =  pu  Ip 

(3.47) 

u"  p'u' 

(3.48) 

u 

pu 

v"  p'v' 

(3.49) 

u 

pu 

pu"v"  fs  pu'v' 

(3.50) 

The  Favre  form  of  the  turbulent  shear  stress  is  related  to 
by  the  following  expression 


cross-correlations  in  velocity 


(3.51) 
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IV.  Results  and  Discussion 


This  chapter  details  the  results  achieved  from  the  mapping  of  the  flowfield 
about  the  wrap-around  fin  (WAF).  A  coordinate  system  is  defined  from  the  nozzle 
throat,  shown  in  Figure  4.1.  All  plots  in  this  section  are  based  on  this  coordinate 
system  and  are  non-dimensionalized  by  the  fin  radius  (r  =  1.5875cm). 

4.1  Inlet  Plane 

The  test  section  in  which  the  WAF  is  tested  is  removable  and  is  mated  to  the 
nozzle  by  a  series  of  alien  screws.  Small  seams  are  present  between  components  of 
the  tunnel,  due  to  its  modular  design.  These  seams  caused  a  series  of  weak  shocks 
and  expansions  which  travelled  into  the  test  section.  Figure  4.2  are  contours  of 
the  mass-flux  in  the  streamwise  direction,  which  shows  a  series  of  these  shocks  and 
expansions  in  the  test  section  inlet.  The  mass-flux  contours  in  this  figure  were  closely 
spaced  about  the  freestream  mean  to  bring  out  the  shock/expansion  detail.  A  sketch 
of  the  inlet  flow  structure  is  shown  in  Figure  4.3.  Although  numerous  shocks  and 
expansions  exist  in  the  inlet,  these  structures  were  weak.  The  negligible  secondary- 
flow  created  by  these  seams  is  shown  in  Figure  4.4  in  the  form  of  mass-flux  vectors. 
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Figure  4.2  Inlet  Plane  Mass-Flux  Contours 


Figure  4.3  Inlet  Seam  Shock/Expansion  Pattern 
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Figure  4.4  Secondary-Flow  at  Inlet 

Note  the  mass-flux  contours  of  Figure  4.2  are,  in  general,  less  than  1.0  when 
non-dimensionalized  by  the  freestream  reference  mass-flux.  The  freestream  values 
are  deflned  based  on  plenum  (settling  chamber)  conditions  and  the  mean  Mach 
number  at  the  test  section  inlet.  The  variation  in  Mach  number  at  the  inlet  to 
the  test  section  is  given  in  Figure  4.5.  The  mean  Mach  number  is  deflned  as  the 
average  of  the  Mach  values  obtained  from  the  center  quarter  of  the  inlet  plane  (p  = 
19.6).  Then  assuming  the  total  temperature  measured  in  the  plenum  to  be  the 
same  value  at  the  test  section  inlet,  the  freestream  static  temperature  is  determined. 
The  static  pressure  is  calculated  using  the  freestream  Mach  number  and  the  Pitot 
pressure.  Using  this  information  in  concert  with  the  equation  of  state  for  a  perfect 
gas  provides  the  freestream  density.  The  freestream  speed  of  sound  is  determined 
from  the  static  temperature,  which  in  turn  yields  the  freestream  flow  speed.  The 
freestream  Reynolds  number  is  then  calculated.  Table  4.1  summarizes  the  freestream 
and  plenum  conditions. 
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Figure  4.5  Inlet  Mach  Contours 


Table  4.1  Tunnel  Inflow  Conditions 


Condition 

Value 

Tt , 

^vlenum 

297K 

Pt  1 

''■plenum 

2.07atm 

T. 

inlet 

113K 

P  1 . 

inlet 

0.065atm 

P  inlet 

0.220^ 

608Si 

s 

Moo 

2.85 

Re 

1.70^ 

_ ID _ 
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4-2  Blending  Region 

A  truncated  section  of  the  model  was  designed  based  upon  WAF  missiles  that 
are  currently  undergoing  flight  tests  at  the  instrumented  test  range  at  Eglin  Air 
Force  Base.  While  a  smaller  scale  missile  of  full  configuration  could  have  been  tested 
in  the  AFIT  supersonic  test  facility,  the  size  of  the  smaller  scale  missile  fins  would 
have  reduced  the  experimental  mapping  resolution.  To  simulate  the  full  missile 
body,  the  aft  section  of  the  missile  was  blended  to  the  tunnel  ceiling.  The  length 
of  the  full  missile  cylinder  was  reduced  by  introducing  a  tapered  region  -  designed 
to  produce  a  weak  shock  coalescence.  The  coalescing  shock  is  visible  in  the  center 
of  the  shadowgraph  of  Figure  4.6  as  the  shock  moving  away  from  the  WAF  model 
(shown  on  the  bottom  of  the  photo). 


Figure  4.6  Composite  Shadowgraph  Image  of  WAF  Model 


Notice  how  the  coalescing  shock  reflected  off  of  the  tunnel  floor  (shown  on  the 
top  of  the  photo)  in  the  schlieren  image  in  Figure  4.7.  Two-dimensional  oblique 
shock  theory  was  used  to  place  the  shock  in  this  position.  The  location  of  the  shock 
was  downstream  of  oblique  shock  theory  and  the  intensity  of  the  shock  was  much 
less  than  theory  due  to  a  smooth  blending  region  and  three-dimensional  relieving 
effects.  The  lower  intensity  shock  results  in  minimal  changes  to  the  flow  qualities. 

Figure  4.8  shows  the  density  field  over  a  cutting  plane  in  the  center  of  the 
blending  region.  The  lighter  region  spanning  the  test  section  was  the  coalescing 
shock  formed  from  the  blending  forebody. 
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Figure  4.7  Composite  Schlieren  Photograph  of  WAF  Model 
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Figure  4.8  Density  Contours  of  Blending  Region 
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Reversing  the  knife  edge  in  a  schlieren  image  produces  the  negative  of  the  pre¬ 
vious  image  (Figure  4.7),  as  shown  in  Figure  4.9.  This  figure  highlights  the  boundary 
layer  on  the  tunnel  floor  (shown  on  the  top  of  the  photo).  A  series  of  shocks  are  also 
visible  extending  from  the  top  of  the  photo  downward.  These  weak  disturbances 
are  created  by  seams  associated  with  removable  plugs.  In  this  schlieren  image  these 
seam  shocks  appear  to  be  much  stronger  than  the  blending  shock.  However,  mea¬ 
surements  presented  in  Figure  4.10  indicate  that  the  plug-seam  shocks  are  of  lower 
strength  than  the  blending  shock  and  are  weaker  than  the  bow  shock.  The  relative 
intensity  of  the  shocks  is  best  visualized  in  the  shadowgraph  photo.  Clearly  the  crisp 
and  intense  character  of  the  bow  shock  suggests  it  provides  the  strongest  influence 
on  the  WAF  flowfield.  The  shock  structure  in  the  test  section  is  summarized  in 
Figure  4.11. 


Figure  4.9  Reversed  Knife-Edge  Schlieren  of  WAF  Model 


4-3  Wrap-Around  Fin  Region 

A  majority  of  the  data  collection  process  for  this  experiment  was  involved  in 
mapping  the  flowfield  in  close  vicinity  to  the  WAF.  The  hot-film  probes  used  for  this 
effort  were  very  fragile.  The  probe  was  destroyed  if  it  brushed  against  the  model. 
When  the  tunnel  was  in  operation,  the  hot-film  probe  would  flex  more  than  1mm 
vertically  (towards  the  model).  To  avoid  probe  contact,  the  probes  were  initially 
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Figure  4.10 


Comparison  of  Shock  Intensity  with  Spanwise  Mass-Flux  Contours 
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Figure  4.11  Overall  Shock  Structure 


positioned  at  a  minimum  distance  of  2mm  from  the  model.  This  limitation  made 
it  impossible  to  examine  the  boundary  layer  near  the  WAF  or  the  juncture  of  the 
WAF  to  the  half-cylinder  body. 

4-3.1  Fin-Body  Juncture.  Despite  the  lack  of  quantitative  measurements  at 
the  intersection  of  the  fin  to  the  forebody,  a  variety  of  qualitative  measurements  were 
made.  A  detailed  shadowgraph  image  is  shown  in  Figure  4.12.  The  A-shock  structure 
visualized  in  the  shadowgraph  has  been  seen  in  similar  juncture  studies  (8,  77,  29). 
Figure  4.13  is  a  schematic  of  the  flowfield  at  the  fin-body  juncture. 


Figure  4.12  Magnified  View  of  Shadowgraph  Near  the  WAF 
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Figure  4.13  Shock  Structure  at  Fin-Body  Juncture 


Fomison  (29)  and  Williams,  HarlofF  and  Gessner  (77)  looked  at  the  vortex 
structure  about  fin-like  junctures  with  flat  plates  as  illustrated  in  Figures  4.14  and 
4.15  from  Fomison.  Oil  flow  visualization  techniques  suggest  that  one  vortex  exists 
in  the  separation  region  under  the  A-shock  on  the  missile  body.  Figures  4.16  and 
4.17  are  photos  taken  of  the  oil-flow  streakline  visualization.  Patterns  seen  in  these 
photographs  show  the  vortex  on  the  convex  side  of  the  WAF  being  much  farther 
away  from  the  model  centerline  than  the  vortex  on  the  concave  side.  The  postulated 
vortex  structure  is  diagramed  in  Figure  4.18. 

4-3.2  Bow  Shock  Region.  The  flowfield  surrounding  a  WAF  fin  differs 
significantly  from  the  flow  structure  present  about  straight  fins.  An  asymmetric  bow 
shock  envelops  the  WAF  as  shown  in  Figure  4.19.  The  structure  on  the  concave 
side  of  the  fin  is  markedly  different  than  the  flow  structure  on  the  convex  side. 
Figures  4.20  and  4.21  show  this  structure  in  Mach  number  contours  progressing  past 
the  fin.  The  convex  side  of  the  fin  is  marked  by  a  large,  sweeping  shock  followed  by  an 
expansion.  A  series  of  contour  composite-images  provide  a  better  overall  description 
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Figure  4.14  Blunt  Fin,  Fin-Juncture  Flow  (Fomison,1986) 
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Figure  4.16  Oil  Flow  Visualization  on  the  WAF  Model,  Concave  Side 
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Figure  4.17  Oil  Flow  Visualization  on  the  WAF  Model,  Convex  Side 


Figure  4.18  Cartoon  of  the  Vortex  Structure  about  the  Fin 
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Figure  4.19  Bow  Shock  Structure 


of  the  shock  structure,  shown  in  Figures  4.22,  4.23  and  4.24.  Figure  4.22  shows  how 
the  Mach  number  drops  on  the  compressive  side  of  the  fin,  while  the  flow  accelerates 
through  the  expansion  in  Figure  4.23.  The  bow-shock  is  shown  from  multiple  cuts 
of  mass-flux  contours  in  Figure  4.24.  Finally,  vertical  cutting  planes  of  streamlines 
show  the  shock  progression  from  the  tip  toward  the  base  in  Figure  4.25. 

A  numerical  study  by  Abate  and  Cook  (9)  proposed  that  the  convex  side  of 
the  WAF  expands  the  flow,  while  the  concave  side  compresses  the  flow.  The  current 
experimental  effort  supports  that  conclusion.  Figure  4.26  contours  the  static  pressure 
in  a  progression  of  planes,  moving  downstream.  The  expansion  shows  up  in  subplot 
(d)  vividly  as  the  low  pressure  section  running  the  length  of  the  convex  side  of  the 
fin.  The  compression  of  the  concave  side  is  illustrated  in  all  of  the  subfigures  as  the 
elliptical  high  pressure  region  just  outboard  of  the  fin. 

The  shock  on  the  concave  side  of  the  WAF  stays  very  close  to  the  WAF  surface, 
while  the  shock  on  the  convex  side  moves  out  of  the  probe  mapping  range  quickly. 
Numerical  studies  of  similar  WAF  configurations  show  good  agreement  with  this 
structure  (7,  9,  70).  In  particular,  Tilmann  (70)  contains  an  inviscid  numerical  sim¬ 
ulation  of  this  WAF  model  run  under  similar  freestream  conditions.  Figure  4.27  is  a 
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Figure  4.20  Mach  Number  Contours 
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Figure  4.21  Mach  Number  Progression 
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r/r=2.6 


Figure  4.22  Concave  Mach  Contour  Map 


Figure  4.24  Mass-Flux  Contour  Map 
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side-by-side  comparison  of  static  pressure  contours,  showing  identical  flow  structure. 
The  agreement  between  the  turbulent,  viscous  experiment  of  this  thesis  to  the  in- 
viscid  numerical  simulation  outside  the  boundary  layer  provides  convincing  evidence 
that  the  bow-shock  structure  is  an  inviscid  phenomenon. 

4.S.3  Tip  Region.  The  bow  shock  which  forms  on  the  beveled  face  of  the 
fin  curves  abruptly  and  travels  over  the  tip  of  the  fin.  The  tip  region  is  marked  by  an 
expansion  that  forms  to  realign  the  flow  interior  to  the  bow  shock.  Figure  4.26  shows 
the  expansion  region  outboard  to  the  WAF  tip,  where  the  static  pressure  is  relieved 
through  the  shock.  Moran  observed  a  tip  bleeding  effect  in  numerical  analysis  of  the 
Hypersonic  Applied  Research  Technology  (HART)  missile  (54).  The  secondary-flow 
streamlines  in  Figure  4.28  do  not  show  similar  trends.  No  bleeding  effect  is  evident  in 
the  secondary-flow  mass-flux  vector  plots  presented  in  Figure  4.29.  Bleeding  effects 
are  also  absent  from  the  measured  flow  angles  about  the  y  and  z  axes,  shown  in 
Figures  4.30  and  4.31.  However,  it  is  important  to  note  the  research  of  Moran  was 
conducted  on  a  missile  model  with  two  straight  fins  at  an  angle-of-attack  (symmetry 
conditions  applied).  Conversely,  the  WAF  experimental  research  conducted  is  on 
a  single  fin  with  no  angle-of-attack.  It  is  expected  that  a  bleeding  effect  would  be 
achieved  by  a  WAF  at  an  angle-of-attack. 

4.4  Turbulence  Measurements 

Even  though  the  present  study  was  unable  to  resolve  the  turbulent  boundary 
layers  about  the  fin,  freestream  turbulence  information  was  obtained.  Presented  in 
this  section  are  the  directly  measured  compressible  turbulent  kinetic  energy  (TKE) 
and  an  estimate  to  the  conventional  incompressible  TKE.  These  data  may  prove  use¬ 
ful  for  CFD.  Especially  important  are  the  effects  of  the  “inviscid”  flow  phenomenon 
on  the  local  freestream  TKE.  These  effects  are  important  because  the  freestream 
turbulence  levels  significantly  influence  the  wall  boundary  layer  and  vis- versa.  Eig- 
ure  4.32  presents  the  incompressible  and  compressible  turbulent  kinetic  energy.  The 
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Figure  4.26  Static  Pressure  Contours 
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Figure  4.27  Comparison  with  Inviscid  Numerical  Solution  (AIAA  96-0190) 


contour  plot  on  the  left  is  the  incompressible  TKE,  while  compressible  TKE  is  on 
the  right.  The  TKE  plots  are  almost  identical  in  the  region  near  the  bow  shock  and 
near  solid  boundaries. 

The  information  in  this  chapter  is  a  summary  of  the  information  included 
in  Reference  (38).  This  reference  provides  a  detailed  summary  for  every  variable 
available,  presented  as  cutting  plane  contours. 
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Figure  4.28  Secondary-Flow  Streamlines  at  WAF  Tip 


Figure  4.29  Secondary-Flow  Mass-Flux  Vectors  at  WAF  Tip 
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Figure  4.30  Flow  Angle  About  Y-Axis 


Figure  4.31  Flow  Angle  About  Z-Axis 
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Figure  4.32  Turbulent  Kinetic  Energy  Comparisons 


V.  Conclusions  and  Recommendations 


The  following  section  summarizes  the  significant  conclusions  following  the  anal¬ 
ysis  of  the  flowfield  in  the  vicinity  of  a  wrap-around  fin  (WAF). 

5.1  Conclusions 

1.  A  high  quality  dataset  was  created  for  numerical  validation.  All  of  the  ap¬ 
plicable  criteria  set  forth  by  Settles  and  Dodson  (61)  have  been  met.  The 
tunnel  inlet  conditions  are  documented  for  use  as  numerical  boundary  condi¬ 
tions.  Mean  flow  and  turbulent  fluctuations  are  presented  in  a  high-resolution 
mesh  in  the  vicinity  of  the  wrap-around  fin.  These  mean  and  turbulent  mea¬ 
surements  are  outside  of  the  fin  boundary  layer,  therefore  the  measurements 
will  not  validate  numerical  turbulence  models.  However,  the  measurements  are 
excellent  for  inviscid  numerical  comparisons  and  freestream  turbulence  infor¬ 
mation.  The  results  also  include  a  detailed  uncertainty  analysis  to  provide  a 
reasonable  bound  on  the  information  provided  in  this  effort. 

2.  The  bow-shock  structure  present  in  the  vicinity  on  a  WAF  is  an  inviscid  phe¬ 
nomenon.  The  experimental  results  of  this  work  were  compared  to  a  compan¬ 
ion  study,  wherein  the  Euler  equations  were  solved  numerically.  The  measured 
and  computed  pressure  contours  were  in  excellent  agreement.  This  suggests  the 
Euler  calculations  may  be  sufficient  to  capture  the  shock  structure  associated 
with  a  full  WAF  missile. 

3.  The  bow-shock  is  asymmetric  about  the  WAF  with  stronger  shock  structure 
on  the  concave  side.  Mean  flow  measurements  were  analyzed  for  differences 
between  the  concave  and  convex  sides  of  the  WAF.  Mach  number  and  static 
pressure  contours  show  that  the  bow-shock  is  much  stronger  on  the  concave 
side  than  the  convex  side.  A  compression  of  the  flow  occurred  on  the  concave 
side  of  the  fin,  while  the  flow  was  expanded  on  the  convex  side.  From  oil-flow 
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streakline  visualization,  it  is  evident  that  the  bow-shock  is  located  farther  away 
from  the  convex  side  of  the  fin  than  from  the  concave  side. 

4.  The  bow-shock  dissipates  over  the  fin  tip.  Through  comparison  of  Mach  number 
and  pressure  contours,  the  shock  disappears  outboard  of  the  WAF  tip.  The 
structure  in  this  region  is  dominated  by  the  expansion  coming  from  the  tip  and 
realigning  the  flow  downstream  of  the  initial  bow-shock. 

5.  A  X-shock  forms  at  the  fin-body  juncture.  Similar  supersonic  studies  of  fins  and 
struts  show  a  shock  “foot”  which  is  pushed  upstream  of  the  fin-body  juncture 
in  the  body’s  boundary-layer.  The  result  of  this  structure  is  a  coalescing  series 
of  shocks,  shaped  like  the  Greek  symbol  A. 

6.  A  single  vortex  structure  at  the  fin  base  is  proposed.  Visualization  of  the  surface 
flow  revealed  a  separation  line  and  a  re-attachment  line  at  the  body,  near  the 
WAF.  The  structure  in  the  separated  region  was  postulated  to  be  consisting 
of  a  single  vortex  which  stays  near  the  body. 

5.2  Recommendations 

The  author  believes  the  current  effort  to  be  only  the  first  step  in  a  series  of 
experiments  aimed  at  characterizing  supersonic  and  hypersonic  curved  body  flow. 
The  current  work  was  based  on  a  DoD  need  to  understand  the  flow  structure  about 
curved  geometries  applied  to  WAFs.  However,  the  general  issue  is  much  larger  than 
missile  applications.  As  construction  techniques  advance,  the  mating  of  curved  sur¬ 
faces  with  other  curved  surfaces  will  become  more  aggressive  on  high  speed  vehicles 
for  the  benefits  of  stealthiness,  aerodynamics  and  stowability.  To  reach  these  goals, 
the  author  suggests  the  following  as  “next-steps.” 

1.  Continue  further  study  on  the  current  WAF  model.  This  study  was  only  the 
first  look  at  a  geometry  of  this  type.  Further  research  needs  to  be  accomplished 
in  definition  of  the  structures  close  to  the  WAF  and  the  fin-body  juncture.  To 
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this  end,  the  model  could  be  mounted  on  the  side  wall  of  the  tunnel,  as  shown  in 
Figure  5.1.  A  side-mounting  of  the  model  is  the  easiest  way  to  reach  the  recesses 
of  the  WAF  with  the  AFIT  Mach  2.9  test  facility.  The  arrows  shown  in  the 
figure  show  the  necessary  probe  locations  which  will  produce  better  definition 
of  the  boundary  layer  and  fin  recesses.  In  addition  to  investigating  model  in 
this  arrangement,  the  boundary  layer  on  the  side  walls  must  be  examined. 
The  effect  of  a  different  boundary  layer  development  could  limit  the  utility  of 
comparisons  to  the  current  effort. 


Figure  5.1  Side-Mounted  WAF  Model  and  Probe  Stations  of  Interest 

Another  research  effort  should  be  conducted  using  a  normal-wire.  Figure  5.2 
identifies  probe  stations  which  will  add  value  to  the  current  effort.  The  normal- 
wire  probe  can  get  much  closer  to  the  solid  boundaries  than  a  cross- wire  probe. 
This  probe  can  better  define  the  boundary  layer  on  the  body,  just  before  the 
juncture  with  the  fin,  and  look  at  the  A-shock  in  detail. 

Additionally,  measurements  of  surface  pressure  should  be  made  on  the  WAF. 
Static  pressure  taps  can  be  drilled  at  positions  from  10%  to  100%  of  the  fin’s 
chord  length  along  the  entire  span  of  the  WAF.  Two  extra  fins  will  be  required 
since  the  fin  is  not  thick  enough  to  run  tubing  through  the  center  of  the  WAF. 
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The  flexible  tubing  will  have  to  be  recessed  into  the  side  opposite  the  pressure 
taps.  Disturbances  created  from  the  tubing  will  limit  the  taps  to  one  side  of 
the  fin. 


Figure  5.2  Normal  Wire  Probe  Stations 

2.  Study  variations  of  this  model.  A  baseline  case  for  the  blending  body  and 
cylinder  is  required.  Then  an  intermediate  step  between  the  baseline  and  the 
current  effort  suggests  a  look  at  the  current  body  with  a  straight  fin.  These 
exercises  will  add  to  numerical  validation  by  providing  steps  for  a  numerical 
analyst  to  follow  while  building  their  model. 

After  baselining  this  effort,  extended  research  on  the  WAF  should  be  aimed 
at  eliminating  the  author’s  simplifications.  The  current  work  looked  at  only 
one  fin,  with  no  angle-of-attack  or  cant  angle,  on  a  partial  body  with  single 
overheat  reduction  techniques.  The  first  effort  might  eliminate  the  question 
of  total-temperature  fluctuations  by  providing  a  cursory  look  using  multiple- 
overheat  normal-wire  reduction.  Multiple  fin  effects  would  best  be  studied 
using  a  semi-cylindrical  model  and  two  WAFs.  Another  variation  would  put 
a  fin  at  varying  angles-of-attack.  Finally,  comparisons  need  to  be  made  with 
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full  missile  mock-ups  using  four  WAFs  to  ensure  the  simplifications  did  not 
eliminate  any  structures. 

3.  Apply  numerical  simulation.  One  of  the  goals  of  this  experimental  research  is 
to  provide  field  data  which  can  be  compared  to  numerical  developments.  Then 
the  envelope  of  information  can  be  drastically  expanded  by  perturbations  in 
Mach  number,  angle-of-attack,  fin  thickness,  fin  sweep  and  orientation  and 
missile  spin. 
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Appendix  A.  Error  Analysis 


This  section  focuses  on  the  known  errors  in  all  measurements  during  the  ex¬ 
periment.  The  effect  of  measurement  errors  is  then  propagated  through  the  data 
reduction  process.  Equations  used  in  data  reduction  were  linearized  to  provide  ap¬ 
proximate  error  bounds  on  the  processed  data  for  this  experiment. 

The  I2  norm  is  utilized  to  combine  errors.  The  I2  norm  is  defined  by 

e-  =  I|e.il2  = 

The  Li  norm  (a  summation  of  the  absolute  error  values)  has  been  found  to  be 
too  conservative  an  estimate  when  compared  to  a  perturbation  analysis  of  the  data 
reduction  equations. 

A.l  Measurement  Errors 

Every  measurement  made  has  an  associated  error  with  it.  It  is  assumed  all 
measurement  errors  are  random,  with  a  Gaussian  distribution.  The  assumption  of 
random  errors  precludes  the  existence  of  biased  errors  or  blunders. 

Pressure  was  measured  by  pressure  transducers,  which  measure  gauge  pres¬ 
sure.  An  error  of  ±0.0017atm  was  accrued  in  measurement  of  ambient  pressure, 
determined  to  be  one-half  of  the  smallest  gauge  increment.  Additionally,  the  man¬ 
ufacturer  of  the  pressure  transducers  has  advertised  an  error  for  Pitot,  cone-static 
and  plenum  pressure  to  be  0.4%,  0.4%  and  0.5%,  respectively  (5).  Errors  in  pressure 
measurements  due  to  calibration  were  determined  to  be  one  standard  deviation  off 
their  respective  calibration  curves.  The  standard  deviation  of  the  calibrations  are 
0.2%  (plenum)  and  0.9%  (Pitot,  cone-static).  Finally,  digital  conversion  of  the  pres¬ 
sure  data  for  storage  and  processing  adds  digitization  error.  All  data  was  digitized 
with  a  12  bit  word  length  (4).  Pressure  measurements  were  sampled  with  a  12.0V 
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range,  resulting  in  an  error  of  ±0.003V.  Multiplication  of  this  voltage  error  by  the 
slope  of  each  calibration  produces  an  error  in  pressure  measurement.  This  error  cor¬ 
responds  to  an  error  of  ±0.0041atm  in  plenum  pressure  and  ±0.00046atm  in  Pitot 
and  cone-static  pressure. 

The  temperature  in  the  plenum  was  measured  during  a  run  at  the  beginning 
of  a  test  period.  This  total  temperature  was  found  to  vary  no  more  than  ±2.36K 
during  subsequent  tunnel  runs  in  that  test  period.  Also,  by  multiple  overheat  (MOH) 
hot-wire  anemometry  in  this  test  facility,  McCann  (51)  and  Miller  (52)  both  found 
boundary  layer  total  temperature  root  mean  square  (RMS)  fluctuations  of  2.0%.  The 
RMS  fluctuations  are  included  as  errors  in  this  analysis  as  they  provide  a  definitive 
bound  of  measurement  uncertainty  from  large  data  samples. 

The  location  of  the  various  probes  in  the  tunnel  was  determined  by  aligning 
two  grids  onto  the  Plexiglas  plates  to  measure  x  and  y  position.  This  measurement 
determined  flexing  of  each  type  of  probe  during  tunnel  operation.  The  flex  angle  and 
position  were  then  determined  as  a  function  of  measured  position.  The  measurement 
of  x  and  y  position  was  accurate  to  ±0. 00127m. 

Without  Plexiglas  windows  on  the  top  and  bottom  of  the  tunnel,  it  was  more 
difficult  to  determine  the  z  position  of  each  probe.  The  error  in  z  position  was 
assumed  to  be  twice  the  error  in  x  or  y.  An  additional  error  from  digitization 
occurred  in  the  measurement  of  y  position  -  as  y  position  was  recorded  using  a 
linear  displacement  voltage  transducer,  LDVT.  Since  the  LDVT  was  sampled  with  a 
12.0V  range,  the  above  error  in  digitization  voltage  was  multiplied  by  the  calibration 
slope.  This  manipulation  resulted  in  an  error  in  y  position  from  digitization  of 
±0. 000015m.  Errors  in  position  were  compounded  by  the  flexing  of  the  probes  during 
tunnel  operation.  Though  the  flex  effect  was  accounted  for,  errors  in  y  measurement 
produced  an  additional  error  in  x. 

Hot-wire  measurements  have  shown  1.0%  root  mean  square  fluctuation  in  volt¬ 
age  for  freestream  flow.  Additionally,  a  digitization  error  of  ±0.003V  was  accounted 
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for  in  the  12.0V  sampling  range.  The  errors  in  pressure,  temperature,  position  and 
voltage  have  been  summarized  in  Table  A.l. 

Table  A.l  Measured  Error  Bounds 


Normalized  Error 

Value 

0.4% 

4.1% 

8.0% 

2.1% 

^Pn 

0.5% 

^Pf>. 

1.1% 

^Pr, 

2.1% 

_ Lin _ 

1.0% 

A. 2  Error  Propagation 

The  errors  listed  in  Table  A.l  have  an  influence  on  all  subsequent  data  reduc¬ 
tion.  To  determine  the  influence  on  calculations,  the  equations  used  are  linearized 
about  freestream  conditions.  Table  A. 2  lists  the  freestream  conditions  for  analysis. 
These  freestream  conditions  were  also  used  to  normalize  the  results. 


A. 2.1  Mean  Flow.  The  majority  of  mean  flow  calculations  are  based  on  the 
local  Mach  number.  The  Mach  number  is  calculated  from  a  curve  fit  to  experimental 
data.  The  curve  fit  is  based  on  the  ratio  of  pressures  from  a  Pitot  probe  and  a  10° 
cone  static  probe.  The  variable  ^  is  introduced  as  the  curve  fit  variable  and  takes 
the  following  form 


p 


Ptu 


P  P  (A.2) 

Ttlcs  -rt2 

The  plenum  pressures  (Ptics,  Ttit2)  were  introduced  to  minimize  errors  caused  by  the 
differences  between  the  two  separate  tunnel  runs,  one  for  each  probe.  Since  ^  was  a 
combination  of  pressure  measurements,  the  errors  combine  with  the  I2  norm 


’  ^Pt2  ’  ^Pc. 


(A.3) 
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Table  A. 2  Freestream  Conditions 


The  curve  fit  (Equation  3.4)  is  linearized  and  provides  an  error  in  Mach  number  of 

=  \\M(C^(  +  2C2e  +  SCsf  +  ^  0.735c^  (A.4) 

The  Mach  number  and  plenum  temperature  were  used  to  derive  the  local  tunnel 
temperature,  Ti.  The  underlying  assumption  to  this  analysis  was  isentropic  flow. 
Although  isentropic  flow  was  violated  as  shocks  and  turbulence  were  encountered, 
no  method  for  measuring  local  temperature  was  readily  available.  The  additional 
assumption  of  an  ideal  gas  and  the  equation  of  state  allowed  for  the  determination  of 
density  and  speed.  Table  A. 3  has  summarized  the  propagation  errors  in  mean  flow 
data  analysis. 


Table  A. 3  Mean  Flow  Error  Bounds 


Normalized  Error 

Derivation 

Value 

2 

2.4% 

0.735e, 

1.8% 

2 

3.1% 

Cq 

O.Sey 

_ _ _ _ _ 

1.6% 

2 

3.6% 

M  -^1 

2 

4.8% 

Cu 

2 _ 

2.4% 

Cy 

2.4% 

2.4% 

A. 2. 2  Turbulent  Flow.  Single  overheat  (SOH)  hot-film  anemometry  was 
used  to  measure  mass-flux  mean  flow  and  RMS  fluctuations  (SOH  analysis  assumes 
negligible  total  temperature  fluctuations  and  has  proved  to  be  valid  for  this  experi¬ 
mental  facility  (52,  24).  Additionally,  cross- wire  measurements  provided  directional 
information.  A  detailed  explanation  of  hot-film  data  reduction  techniques  is  pre¬ 
sented  in  Chapter  III.  Due  to  the  complicated  structure  of  hot-film  data  reduction, 
a  logarithmic/derivative  technique  was  applied  to  the  hot-film  data  reduction  equa¬ 
tions  to  estimate  the  propagation  of  errors  throughout  hot-film  analysis. 
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The  Nusselt  number  (a  measure  of  heat  transfer  from  an  infinite  cylinder),  Nu, 
was  determined  from  the  power  consumption  required  to  maintain  a  constant  wire 
temperature  (Equation  3.6).  The  error  in  Nu  was  a  function  of  the  measured  voltage 
error  and  total  temperature  error. 

The  effective  Reynolds  number,  Reg,  was  determined  by  a  curve  fit  of  Nu 
(Kings  Law,  Equation  3.7).  Additionally,  the  fluctuation  in  Reg  was  a  function  of 
RMS  voltage.  These  errors  were  combined  under  the  I2  norm. 


Table  A. 4  Turbulent  Flow  Error  Bounds 
Normalized  Error  Derivation  Value 


A. 2. 3  Separation  of  Variables.  In  separating  the  primitive  variables  from 
hot-wire  data,  the  error  in  Reynolds  number  equates  to  an  error  in  mass-flux.  It 
was  assumed  that  errors  in  off-axial  velocities  (ie.  v  and  w)  were  the  same  as  the 
associated  error  in  u.  The  separation  equations  were  linearized  about  the  refer¬ 
ence  conditions  of  Table  A. 2.  The  results  of  variable  separation  error  analysis  are 
summarized  in  Table  A. 5. 
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Table  A. 5  Variable  Separation  Error  Bounds 


Derivation  Value 


^Rez 

8.1% 

16.3% 

10.3% 


Appendix  B.  Nozzle  Coordinates 


Tables  B.l  and  B.2  detail  the  coordinates  of  the  nozzle  of  AFIT’s  Mach  2.9 
test  facility.  Note  the  nozzle  is  two  dimensional  and  has  a  width  (z)  of  6.35cm.  The 
origin  is  located  on  the  floor  at  the  throat.  The  coordinates  extend  from  upstream  of 
the  convergent  section  to  the  beginning  of  the  first  removable  test  section.  Figure  B.l 
is  a  plot  of  the  data. 
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Table  B.l  Nozzle  Coordinates  (Top) 


y  (cm) 

X  (cm) 

y  (cm) 

X  (cm) 

y  (cm) 

-8.102600 

8.359140 

-1.663700 

3.776980 

1.452880 

1.772920 

-4.292600 

8.359140 

-1.577340 

3.578860 

1.579880 

1.826260 

-4.206240 

8.348980 

-1.488440 

3.385820 

1.645920 

1.855216 

-4.117340 

8.323580 

-1.402080 

3.197860 

2.275840 

2.151126 

-4.030980 

8.282940 

-1.313180 

3.014980 

3.345180 

2.621280 

-3.942080 

8.227060 

-1.226820 

2.839720 

3.848100 

2.824480 

-3.855720 

8.158480 

-1.137920 

2.672080 

4.348480 

3.014980 

-3.766820 

8.074660 

-1.051560 

2.511806 

4.851400 

3.200400 

-3.680460 

7.975600 

-0.962660 

2.360422 

5.367020 

3.380740 

-3.591560 

7.866380 

2.218690 

5.897880 

3.558540 

-3.505200 

7.747000 

-0.787400 

2.086864 

6.449060 

3.733800 

-3.416300 

7.614920 

-0.701040 

1.966214 

7.025640 

3.909060 

-3.329940 

7.472680 

-0.612140 

1.856994 

7.630160 

4.081780 

-3.241040 

7.322820 

-0.525780 

1.759966 

8.265160 

4.251960 

-3.154680 

7.162800 

-0.436880 

1.676146 

8.933180 

4.422140 

-3.065780 

6.992620 

-0.350520 

1.605788 

9.636760 

4.589780 

-2.979420 

6.817360 

-0.261620 

1.549908 

10.38098 

4.754880 

-2.890520 

6.634480 

-0.175260 

1.509268 

11.17092 

4.917440 

-2.804160 

6.446520 

-0.086360 

1.484122 

12.00658 

5.080000 

-2.715260 

6.253480 

1.475740 

12.89304 

5.237480 

-2.628900 

6.055360 

■ilfcWrfdil 

1.482344 

13.83538 

5.387340 

-2.540000 

5.854700 

0.292100 

1.490472 

14.83614 

5.532120 

-2.453640 

5.648960 

0.403860 

1.501648 

15.90294 

5.674360 

-2.364740 

0.513080 

1.516380 

17.03832 

5.803900 

-2.278380 

5.232400 

0.624840 

1.534414 

18.23720 

5.925820 

-2.189480 

5.021580 

0.739140 

1.556004 

19.53260 

6.037580 

-2.103120 

4.810760 

0.853440 

1.581404 

20.92960 

6.136640 

-2.014220 

4.599940 

0.967740 

1.610868 

22.40280 

6.217920 

-1.927860 

4.391660 

1.087120 

1.644396 

23.97760 

6.283960 

-1.838960 

4.183380 

1.206500 

1.682496 

25.65400 

6.329680 

-1.752600 

3.980180 

1.328420 

1.725168 

27.45740 

6.350000 

B-2 


Table  B.2  Nozzle  Coordinates  (Bottom) 


X  (cm)  y  (cm) 


-8.10260 

-6.883400 

-4.29260 

-4.03098 

-6.751320 

-3.85572 

-6.682740 

-3.76682 

-6.598920 

-3.68046 

-6.499860 

-3.59156 

-6.390640 

-3.50520 

-6.271260 

-3.41630 

-6.139180 

^^1 

-3.15468 

00 

1 

(”060 

-3.06578 

-5.516880 

m 

-2.89052 

-5.15^ 

3740 

-2.80416 

-4.970780 

-2.71526 

-4.777740 

X  (cm)  y  (cm) 


-2.10312  -3.335020 


-2.62890 

-2.54000 

-2.45364 

-2.36474 

-2.27838 

-2.18948 


-4.579620 

-4.378960 

-4.173220 

-3.964940 

-3.756660 

-3.545840 


-1.92786 


-1.83896 


-1.75260 


-1.66370 


-1.57734 


-1.48844 


-1.40208 


-1.31318 


-1.22682 


-1.13792 


-1.05156 


-0.96266 


-0.87630 


-0.78740 

-0.70104 


-0.61214 


-0.52578 


-0.43688 


-0.35052 

-0.26162 

-0.17526 

-0.08636 

0.00000 

27.4574 


-2.915920 


-2.707640 


-2.503678 


-2.301748 


-2.103882 


-1.910334 


-1.722120 


-1.540002 


-1.364488 


-1.196086 


-1.036066 


-0.884682 


-0.742950 


-0.611124 


-0.490474 


-0.381254 


-0.284226 


-0.200355 


-0.130100 

-0.074244 

-0.033452 

-0.008481 

0.000000 

0.000000 
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Abstract 

A  ceiling-mounted  semi-cyliiulrical  model  containing  a  single  wrap-around  fin  (WAF)  was  tested  in  the 
AFIT  Mach  2.9  test  facility.  Flow  visnali/.ation  using  oil-flow  streaklines,  schlieren  images  and  shadowgraph 
photography  revealed  a  A-sliock  at  the  fin-body  juncture  and  the  development  of  an  asymmetric  bow-shock 
about  the  fin.  Quantitative  nieasiirements  were  taken  with  a  10°  cone-static  pressure  probe,  a  Pitot  pressure 
probe  and  two  cross-wire  hot-film  probes  (u-\'  and  u-w  components,  respectively).  Measurements  were  made 
at  cutting-planes  from  the  inlet  of  the  lest  section  to  aft  of  the  model,  with  emphasis  placed  in  the  vicinity 
of  the  WAF.  Results  include  cut  I  ing-plane  profiles  and  contours  of  mean  and  turbulent  fluctuations  of  the 
primitive  and  conserved  flow  variables.  It  was  found  that  the  incompressible  turbulent  fluctuating  quantities 
are  equally  as  descriptive  of  the  flow  structure  in  the  fin’s  vicinity  as  the  compressible  turbulence  fluctuations. 
The  asymmetric  bow-shock  was  found  to  be  an  inviscid  phenomenon  which  was  stronger  on  the  concave  side 
than  the  convex  side  and  of  diminishing  strength  at  the  tip  with  no  bleeding  effects  over  the  tip. 
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